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Abstract 

We study recent data on deep inelastic e+p scattering at HERA to constrain 
the parameters of a Randall- Sundrum-type scenario of quantum gravity with 
a small extra dimension and a non-factorable geometry. 
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Theories with extra dimensions which predict observable consequences 
at the current high energy accelerators have lately attracted a great deal of 
interest. Following the original suggestion by Arkani-Hamed, Dimopoulos 
and Dvali (ADD) [|^, there have been numerous studies in the literature |0] 
which probe consequences of multiple Kaluza-Klein graviton exchange lead- 
ing to interactions of electroweak strength. The fact that these theories 
predict quantum gravity effects at TeV scales has been suggested as a 
solution of the well-known hierarchy problem in the Standard Model (SM). 
Though novel and interesting, however, the model suggested by ADD, which 
is based on a factorable x (S^)'' geometry, d being the number of extra 
compact dimensions, has the drawback of introducing large compactification 
radii (amounting to an energy scale as low as 10"^^ GeV), which effectively 
introduces a new hierarchy problem. Motivated by this, Randall and Sun- 
drum (RS) have suggested a somewhat different mechanism to solve the 
hierarchy problem. Instead of writing a factorable metric 

ds'^ = ri^j,y dx^ dx" + Rl d(j)i d(f)i (1) 

where the 0j {i = l,d) are extra dimensions compactified with a common 
radius ~ 1 mm, they write a non-factorable metric 

ds^ = e-^^^'^ 7]^^ dx^'dx" + Rl d(j)^ (2) 

involving one extra dimension compactified with a radius Rc, which is as- 
sumed to be marginally greater than the Planck length 10^'^^ cm, and an 
extra mass scale /C, which is related to the Planck scale Mp in the five- 



dimensional bulk by /C Mp^ ~ Mp' . Such a 'warped' geometry is 
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motivated by compactifying the extra dimension on a S^/Z2 orbifold, with 
two D-branes at the orbifold fixed points, viz., one at = ('Planck brane' 
or 'invisible brane'), and one at = vr ('TeV brane' or 'visible brane'). It 
can then be shown that if we assume matter fields to be confined to these 
Z^-branes, one can solve the Einstein equations to obtain a metric of the 
above form. The interesting physical consequence of this geometry is that 
any mass scale Ai on either brane gets scaled by the 'warp factor e'^^"'^ on 
either brane. Thus, a mass scale on the Planck brane {(f) = 0) will remain 
unchanged, but any mass scale on the TeV brane will be scaled by a factor 
g-7r/ci?,c jf assume that the Planck scale is the only fundamental mass 
scale in the theory, all masses on the TeV brane will be scaled to 

M ~ e-^'^^^M^,') (3) 

It now requires )CRc ~ 11 — 12 to obtain Ai of the order of the electroweak 
scale, which justifies the name 'TeV brane'. Thus, in this model there is 
no hierarchy problem, since all the independent mass scales are close to 
the Planck scale. There still remains a minor problem: that of stabilizing 
the radius Rc (which is marginally smaller than the Planck scale) against 
quantum fiuctuations, but this is not so severe as in the model of ADD, 
where the compactification radius needs to be stabilized over as many as 30 
orders of magnitude. A simple extension of the RS construction involving an 
extra bulk scalar field has been proposed to stabilize Rc and this predicts 
light radion excitations with possible collider signatures However, as 
these will not contribute to the processes of interest in this letter, this idea 
will not be discussed further. On the fiip side, it is not as simple to embed 
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the RS construction within the framework of string theories as it is for the 
ADD case. However, a first attempt has been made [Q, and it may be hoped 
that future work will achieve this highly desirable goal. 

Following the ingenious suggestion of a non-factorable geometry, the mass 
spectrum and couplings of the graviton in the RS model have been worked 
out, in Refs. 0, |[. We do not describe the details of this calculation, but 
refer the reader to the original literature. It is worth noting that there are 



strong phenomenological constraints on bulk excitations of the SM fields 11 



It suffices here to note that the effective Lagrangian density for graviton 
interactions on the TeV brane (which we identify with the observable world) 
has the form [Bl 



1 pirKRc 

^f// = -jT ^J^(^) ^'^(^) -J^T. h;A^) T^-ix) (4) 

IVlp Mp ^^-^ 



where Mp = Mp/VSn is the reduced Planck mass and the /i" (x) corre- 
spond to the Kaluza-Klein expansion of the massless graviton in five dimen- 
sions 

K.{xA) = E h;^{x)^. (5) 

n=0 V -ttc 

Equation (^) tells us that the massless Kaluza-Klein (KK) mode effectively 
decouples from ordinary matter since its interactions are suppressed by the 
Planck mass. On the other hand, the massive KK modes couple as the 
inverse of the Planck mass, scaled by e~'^^'^=, which is an electroweak-strength 
interaction. Feynman rules to the lowest order for these modes, assuming a 
coupling Mp^ to all the modes have been worked out in Refs. [^]) and |1^ in 
the context of ADD-like scenarios. All we need to do to get the corresponding 
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Feynman rules in the RS model is to multiply the couplings by the warp factor 
^ttJCRc ^ff\^QJ^e necessary. 

As shown in Ref. , the orbifold geometry forces the Fourier coefficients 
to satisfy a Bessel equation, whence it may be shown that they are 
given by a linear combination of the Bessel and Neumann functions of order 
2. The requirement that the first derivative of x^i^P) be continuous at the 
orbifold fixed point = vr then requires Ji(x„) = 0. Using this, the masses 
M„ of the graviton states can be written in terms of the zeros of the Bessel 
function of order unity as 

M„ = x^/Ce-"'^^^ = x^rrio (6) 

where mo sets the scale of graviton masses and is essentially a free parameter 
of the theory. It is also convenient to write 

p-wKRc ^ 

= ^ (7) 

Mp mo ^ ^ 

using (|^) and introducing another undetermined parameter Cq = IC/Mp\ 
Ref. 1^ points out that (mo,Co) may conveniently be taken as the free pa- 
rameters of the theory, and we follow their prescription in our work. 

Though Co and mo are not precisely known, one can make estimates of 
their magnitude using theoretical ideas and phenomenological inputs. We 
note that the RS construction requires /C to be at least an order of magni- 
tude less than Mp \ because /C~^ sets the scale for the curvature of the fifth 
dimension, and should therefore be large compared with the Planck length. 
The latter is necessitated by the requirement that fiuctuations in the bulk 
gravitational field in the vicinity of the D-branes be small. The range of 



interest for cq is, therefore, about 0.01 to 0.1, the lower value being deter- 
mined by naturalness considerations. Regarding mo, Eq. (|) tells us that it 
is reduced from the scale JC by the factor e~'^^^''. In the RS construction, one 
requires )CRc ~ 11-12, which reduces ttiq to the electroweak scale. Hence, 
we may consider mo in the range of a few tens of GeV to a few TeV. Eq. 
also tells us that the first massive graviton lies at Mi = ximo — 3.83 mo- 
Since no graviton resonances have been seen at LEP-2, running at energies 
upto 200 GeV, it is clear that we should expect mo > 52 GeV. 

In this letter we report on a study of graviton effects, within the RS model, 
on e~^p deep inelastic scattering (DIS) at HERA. At the leading order, there 
are two extra Feynman diagrams contributing to e^p — > + X. One of 
these involves a t-channel exchange of a virtual (massive) graviton between 
the e"*" and a quark; the other involves a t-channel exchange of a virtual 
(massive) graviton between the e"*" and a gluon. The first one adds coherently 
with the corresponding SM diagrams with photon and Z-boson exchange; the 
second one has no SM analogue and hence adds incoherently. However, at 
HERA energies, we do not expect much contribution from the gluon-induced 
diagram because of the low gluon flux. 

The cross-section for the above processes has been calculated for the case 



of the ADD model in Ref. |]12[ and can be easily translated to the RS model 
using the replacement 

Ml ml ^\i\+M^ ^ ' 

We have developed approximate analytic formulae for this sum, using the 
well-known properties of the zeros of the Bessel function Ji {x) . These will be 



presented elsewhere |]T3 . Using this, we incorporate the calculated theoretical 



cross-section into a parton-level Monte Carlo event generator, with two free 
parameters, viz. the graviton mass scale mo, and the coupling parameter 
Cq. Finally the simulation results are compared with data from the ZEUS 
Collaboration to constrain the {rriQ-Co) plane. 

Our numerical studies are founded on the latest results presented |T^ by 
the ZEUS Collaboration, which are based on 47.7 pb^^ of data collected over 
the period 1994-1997. The ZEUS Collaboration uses the double-angle (DA) 
method to determine the DIS variables. In this method, one measures the 
polar angle 9e of the scattered positron, and reconstructs the polar angle ■jh 
of the struck quark in the naive parton model using all hadronic clusters 
which can be identified with a jet having the requisite pt balance with the 
positron. In terms of these observables and the energy Ee{Ep) of the initial 
positron (proton) beam, one can reconstruct the standard DIS variables as 

n2 Sin7fe(l + COSge) .Q^ 

Qda = ^^e- \ — ^ —( (yj 

sm + sm Oe - sm(7/, + 0^) 

Ee sin + sin 9^ + sm{jh + de) ^.s 
Ep sm 'jh + sm 9^ - sm{'jh + 9e) 

sm9e{l - cosjh) ,,,s 
sm -fh + sm 9e - sm[-fh + 9e) 

Various triggers, acceptances and selection cuts have been used by the ZEUS 
Collaboration, of which we need to impose only the following in a parton-level 
analysis: 

• If the final state positron has polar angle greater than 17.2^^, it must 
have a total energy greater than 10 GeV. 
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• If the final state positron has polar angle less than 17.2°, it must have 
a transverse momentum greater than 30 GeV. 

Of the DIS variables listed above, it is well-known that it is the first, 
namely Q'oa, which exhibits maximum sensitivity to most kinds of new 
physics. The ZEUS Collaboration has presented their data for 20 bins in 
Qj^A^ ranging from Qj^^ = 400 GeV^ to 51200 GeV^. The Born-level cross- 
section in each bin, obtained by suitably scaling out radiative effects, has been 
presented by the ZEUS Collaboration together with the SM expectation. We 
have checked that the latter, obtained using hadronization procedures incor- 
porated in the standard HERACLES and ARIADNE program packages, are 
in excellent agreement (within a few per cent) with our parton-level analysis. 
Any small differences which persist can be removed by calibrating the cross- 
section binwise, so as to yield the actual ZEUS expectations. This procedure 
has the added merit of taking care of residual higher order effects such as 
initial state radiation, which should be roughly the same in the SM as in the 
case when the graviton exchanges are included. 

In Fig. 1, we present a graph showing the variation in the Q%j^ distri- 
bution with the mass scale mo in the RS model. For this graph, we have 
plotted the ratio 

R{mo,co) = (12) 

of the cross-section predicted in the RS model with the prediction of the SM, 
for Co = 0.1 and rriQ = 80, 90 and 100 GeV, together with the ZEUS data. It 
may be seen that the cross-section in the RS model (like the ADD model |T2| ) 



exhibits large deviations in the highest Q bins. This, of course, rapidly 



approaches the SM (dotted hne) if Cq is chosen smaller. Interestingly, the RS 
model predictions seem to show a slight diminution for intermediate values 
of Qj^A^ which are intriguingly like the trend shown by the data. However, 
the experimental errors are too large to enable us to attach any significance 
to this circumstance. Accordingly, we take the conservative viewpoint that 
the data fit the SM very well and can be used to constrain new physics. 
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Figure 1. Illustrating the ratio i?(mo,co) of the Q%y^ distribution in the RS model to 
that in the SM (see Eq. [7^, for cq = 0.1 and uiq ~ 80,90 and 100 GeV. The dotted line 
corresponds to the SM. The ZEUS data are also shown. 

Once the above simulation is set-up, we estimate the binwise cross-section 
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X^if^o, Co) for each value of tjiq and Cq and use this distribution to calculate 

2 



X^(mo,co) = J2 



20 [(Ti(mo,co) - af'^'^ 



(13) 



i=l 



where 



e, = el ^[a,(mo,co)-af + ^[af''-) - a,(mo, cq)] (14) 

assuming that the experimental value in the z-th bin is given by [cf'"'^'']^^, . In 
this, it is assumed that e\ and e\ contain the statistical and systematic errors 
added in quadrature. The 95% C.L. bound is then obtained by requiring 



X (?Tio;Co) < 31.41, which is the expectation from random fluctuations. 
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Figure 2. Illustrating the constraint on the parameter space of the RS model arising 
from an analysis of ZEUS high-Q^ data. The shaded region is ruled out at the 95% C.L. 
level. 
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In Figure 2, we show the 95% C.L. constraints on the rriQ-Co plane us- 
ing the above technique. Since the effective graviton couphng is quadratic 
in Co we expect the cross-section to rise as cq increases — this is reffected 
in the fact that Figure 2 shows upper bounds on cq. On the other hand, 
the mo dependence of the cross-section is very comphcated, because of the 
summation over the KK states. However, as the figure makes clear, there 
is a sharp drop in the cross-section as mo increases, so that the ZEUS data 
become insensitive to the new physics beyond about rriQ = 120 GeV. This 
corresponds to Mi ^ 460 GeV, a value which is still not accessible to the 
generation of colliders running at present. 

As the above figure and discussion makes clear, HERA data as presented 
by the ZEUS Collaboration provide somewhat modest, but nevertheless in- 
teresting constraints on the parameter space of the RS model of quantum 
gravity. Since gravitons couple to the energy-momentum tensor of the mat- 
ter fields, one may expect considerable improvements in these results at ma- 
chines running at higher energies, such as the LHC, the proposed NLC and 
possible muon colliders. In particular, it would be interesting to see if these 
machines could actually find graviton resonances, which one expects in the 
RS model |^ , but not in the ADD theory. We have performed a preliminary 
study of the RS model in the light of HERA data, and we expect that the 
future will see many more detailed studies of this very interesting scenario. 
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Dilip K. Ghosh, Sudipta Mukherji, Gautam Sengupta, K. Sridhar and Ajit M. 
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